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This paper presents two-fluid simulations of forced magnetic reconnection with finite electron inertia in a 
coUisionless three-dimensional (3D) cube with periodic boundaries in all three directions. Comparisons are 
made to analogous two-dimensional simulations. Reconnection in this system is driven by a spatially localized 
forcing function that is added to the ion momentum equation inside the computational domain. Consistent 
with previous results in similar, but larger forced 2D simulations [B. Sullivan, B. N. Rogers, and M. A. Shay, 
Phys. Plasmas 12, 122312 (2005)], for sufficiently strong forcing the reconnection process is found to become 
Alfvenic in both 2D and 3D, i.e., the inflow velocity scales roughly like some fraction of the Alfven speed based 
on the upstream reconnecting magnetic field, and the system takes on a stable configuration with a dissipation 
region aspect ratio on the order of 0.15. 



I. INTRODUCTION 

Magnetic reconnection allows a magnetized plasma system 
to convert magnetic energy into high speed flows and thermal 
energy, and is the main driving mechanism behind energetic 
phenomena such as solar flares, magnetospheric substorms, 
and tokamak sawteeth. A defining feature of these phenom- 
ena is that they are very bursty: in a substorm, for example, 
a significant fraction of the lobe magnetic flux is reconnected 
over a period of about 10 minutes, while the time between pe- 
riodic substorms is roughly 2.75 hours Similarly, in the 
case of an X-class solar flare, the energy release may occur 
during several minutes or less in active regions that can last 
for weeks before reconnecting |2|. A viable model of recon- 
nection in such systems must therefore be consistent with the 
relatively long, quasi-stable periods in between reconnection 
bursts, as well as the rapidity of the energy release once a burst 
has somehow been triggered. 

Regarding the rapidity of reconnection, recent research 
(see, e.g., Refs. |l3l|4l|5l|6l|71|8l|9l[l0l[II][I2[l3l[Il[I2 
[T6l [13) has suggested that the reconnection process can in- 
deed generate sufficiently fast rates of energy release to be 
consistent with observations provided that various conditions 
are met. For example, studies of spontaneous, unforced re- 
connection in systems with sufficiently narrow current sheets 
(6 di = c/ojpi in the case without a guide magnetic field) 
have found that the reconnection process becomes Alfvenic - 
that is, the inflow velocity of plasma and magnetic flux into 
the reconnection region scales roughly like some fraction of 
the Alfven speed based on the upstream reconnecting compo- 
nent of the magnetic field. The narrowing of the current layer 
down to widths S ^ di serves as a trigger point for the on- 
set of fast reconnection in this case due to the important role 
played by non-ideal effects such as the Hall term |3 |. Without 
the addition of ad-hoc anomalous resistivity terms, simula- 
tions of systems with substantially thicker current sheets do 
not exhibit such fast reconnection behavior: at best the recon- 
nection rates are throttled by the formation of long, narrow 
outflow layers as in the Sweet-Parker model, or in unfavor- 
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able magnetic geometries (e.g. A' < |fT8|) no spontaneous 
reconnection may take place at all. 

The fact that reconnection can occur very rapidly or not at 
all is consistent with the burstyness of the reconnection phe- 
nomena discussed earlier provided one can explain how a sys- 
tem can transition from a quasi-stable phase (characterized, 
for example, by a relatively thick current sheet) to a period 
of fast reconnection (a narrow current sheet). One possible 
mechanism to explain this transition is explored in this paper: 
to externally force a stable, thick current sheet system so that 
the magnetic flux and electric current are compressed in a lo- 
calized region down to the small scales S ^ di at which non- 
ideal effects such as the Hall term become important. This 
compression is generated in our simulations by a spatially lo- 
calized external forcing term that is added to the ion momen- 
tum equation inside a three-dimensional (3D) computational 
domain. This method allows the scaling of the reconnection 
rate in 3D to be studied over a wide dynamic range: by chang- 
ing the level of forcing, our simulations achieve more than 
an order of magnitude variation in the magnetic field just up- 
stream of the reconnection region, and more than two orders 
of magnitude variation in the reconnection rate. The simula- 
tions show that such localized forcing in 3D can indeed in- 
duce fast rates of reconnection similar to those observed in 
2D, spontaneously reconnecting, narrow current sheet sys- 
tems, even in magnetic geometries that would not reconnect at 
all in the absence of forcing. This result - that a system with a 
locally narrow current sheet will exhibit fast reconnection be- 
havior (whether the narrow current sheet is created by external 
forcing or is present initially) - is reasonable but not obvious, 
since the magnetic geometries outside the forced zone in the 
two cases may be very different. 

This finding extends to 3D a similar result that was obtained 
in a larger 2D study [19], as well as an earlier 2D multi-code 
study of forced reconnection known as the Newton Challenge 
ll24l l25l l26l . The simulations in the Newton Challenge were 
initialized with a relatively narrow current sheet that was un- 
stable to spontaneous reconnection but was somewhat thicker 
than the threshold required for fast (Alfvenic) reconnection 
(two ion skin depths compared to roughly one). The current 
sheet was then pinched by a spatially and temporally depen- 
dent inflow velocity imposed at the upstream boundary of the 
simulation domain. The result was that Alfvenic reconnection 
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indeed occurred in all studies that included the Hall term, al- 
beit with reconnection rates ^ 20% — 50% lower than those 
found in studies of spontaneous (unforced) reconnection in 
systems with somewhat narrower current sheets (one ion skin 
depth). The work described here goes beyond the Newton 
Challenge to explore, over a wide range of parameters, the 
scaling of forced magnetic reconnection in a tearing mode sta- 
ble (A' < 0) system with an initially broad current sheet and 
a forcing function that is spatially localized in 3D. 

Several limitations of the present work that are worth not- 
ing. Our simulations explore the dynamics of forced recon- 
nection in a cubical box of 25.6 ion skin depths (the 3D ana- 
log of the GEM Challenge system |3 1) using a simple isother- 
mal two-fluid model |9| and periodic boundary conditions. 
The frozen-in condition is broken at the smallest scales in 
this model by electron inertia and numerical diffusion effects. 
Particle simulations of reconnection, on the other hand, have 
shown that non-gyrotropic pressure tensor effects, which are 
not included here, should play the dominant role in this re- 
gard. Past 2D studies of unforced reconnection (see, for ex- 
ample, Refs. |[3l|4]|5l|6l|7l[8l) have suggested that the recon- 
nection rate is not strongly sensitive to this difference, and it is 
hoped that this insensitivity carries over to the present study. 
In addition, 2D particle simulations in much larger boxes are 
currently being used to study the sensitivity of the reconnec- 
tion process to the boundary conditions and the physics of the 
electron diffusion region |E1 [28l |29l |30l [iD . Due to the size 
limitations imposed by the 3D nature of our study, as well as 
the physical simplicity of our model, these important issues 
lie outside the scope of the present work. 

The organization of this paper is as follows. In section |ll] 
we describe the simulation model, equilibrium profiles, and 
forcing function. In section [III] we present and discuss the 
results of the simulations included in this study. The main 
conclusions are summarized in section HVl 



II. SIMULATION MODEL 

Our simulations are based on the two-fluid equations with 
the addition of a spatially localized forcing term in the ion 
momentum equation. These equations in normalized form are 



Pe = nTe , Pi ^nTi , p = Pi 



(5) 



B' 



B , Ve = V, - J/n , J = V X B (6) 



Here Fy — 



Fy{x, y, z, t) is the forcing function, the form 
of which is discussed in below. For simplicity we assume 
an isothermal equation of state for both electrons and ions 
(qualitatively similar to the adiabatic case), and thus take 
Te and Ti to be constant. The normalizations of Eqns. ([T]i- 
(|6]l are based on constant reference values of the density 
riQ and the reconnecting component of the magnetic field 
Bxo, and are given by (normalized physical units): t — + 

n/riQ, B 'B/B^Q, V^^e -> 
LUcidi = Ba;o/(47rnomi)i/2, Ti,e 

iir/Blo, J ^ J/(noeVAx), and 



Pi. 



'^r,e/VAx, Vax = 
T.^e^-KTlQ/BlQ, p,, 

Fy Fydi / {m.inoVj^) . Our algorithm employs fourth- 
order spatial finite differencing and the time stepping scheme 
is a second-order accurate trapezoidal leapfrog 132] [33l . We 
consider a cubical 3D simulation box with physical dimen- 
sions Lx Lx L — 25.6di x 25.6di x 25.6di (so that L — 25.6 
in normalized units) and periodic boundary conditions im- 
posed at a; = ±L/2,y — ±L/2, and z = ±L/2 . The sim- 
ulation grid is Tij. X Uy X riz = 128 x 256, xl28, yielding 
grid scales of = A^ = 0.2 and A^ = 0.1. The electron 
to ion mass ratio is me/mi ~ 1/25 so that the normalized 
electron skin depth is de = y^me/nii = 0.2. The frozen-in 
law for electrons is broken primarily by the presence of finite 
electron inertia, which is manifested by the rrie terms in the 
definition of B' in Eq. (|6|. Very near the x-points, however, 
the electron inertia terms become weak in quasi-steady condi- 
tions, since both the dt and v • V terms on the left-hand-side 
of Eq. ([TI become small. In these small regions the frozen-in 
law is mainly broken by numerical diffusion. Past studies of 
this effect (e.g. Ref f34|) have found that the rates of recon- 
nection are relatively insensitive to such grid-level diffusion 
and are in reasonable agreement with the rates obtained from 
particle simulations ll35l . 



n{dtVr+Vi-\7V,) = 3 xB-\/p+Fyy (1) 



dtB' = -V X E' (2) 



A. Initial Equilibrium 

We begin with a one dimensional equilibrium. Three di- 
mensional effects will enter via the forcing function described 
in the next section. The normalized magnetic field and density 
profiles in our initial equilibrium are given by: 







B{y) = sin 


'2tt ( LV 








(J X B' - Wpe) 


(3) 


n{y) = 1 + 






2(T, + Te) 



(7) 
(8) 



dtn + V; • Vn = 



(4) 



The boundary conditions are periodic in all three directions. 
As required by these boundary conditions, B is periodic under 
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y ^ y + L. Note that the normaUzation parameter B^o has 
been chosen so that the peak value of this initial field is unity. 
From the given form of n it is apparent that rt = 1 at this 
location, or in physical units, n = uq. The density profile is 
chosen to satisfy the total pressure balance condition, which 
in normalized form is given by 



n{Ti + Te) + l-Bl^ constant 



(9) 



Unless otherwise stated, we take the (constant) total temper- 
ature to be Ttot = Ti + Te = 1.0, or in physical units 
AirnoTtot / B'^Q — 1, so that the plasma f3 based on the re- 
connecting field has a (minimum) value of 2 where B^ ~ I 
and n = 1. The density reaches a maximum value of n = 1.5 
in the center of the equilibrium current sheets (y ~ ±L/4) 
where B^ = 0. Initially the current is carried entirely by 
the electrons, while the ions begin at rest. To prevent en- 
ergy buildup at the grid scale, the simulations include fourth 
order dissipation in the density and momentum equations of 
the form /i4V^ where /i4 — 5.1 • 10~^. To avoid physically 
artificial effects that can arise from exact reflection symme- 
tries of the initial condition, a small amount of random noise 
is added to the magnetic field and ion current at the levels 



\B„ 



10- 



\J^. 



10- 



The linear tearing mode stability parameter is defined as 
A' = [dy^piy 0+) - dyi){y 0-)]/tP{y = 0) where 4j 
is the perturbation in flux the function 118J . In our system it is 
given by ll36l : 



A'(fc,) = - 2ko^ 



fcg 



1 tan 




(10) 



where fco — 2-k/L. Since periodicity along the x direction re- 
quires kx > 2-k/L and thus k'^/kg > 1, one sees that A' < 0. 
Therefore, as noted in the introduction, the system we con- 
sider is stable to tearing modes in the absence of forcing and 
exhibits no spontaneous reconnection. 



B. Forcing Function 

The initial ID equilibrium described in the previous sec- 
tion does not contain a pre-seeded x-line. Rather, we em- 
ploy a forcing function to drive plasma and magnetic flux into 
the region surrounding the point {x, y, z) ~ {L/4, —L/4, 0), 
thereby forming an x-line along z through that location. 
The forcing function has the general form Fy{x,y,t) = 
X{x)Y{y)Z{z)Q{t). Along x and z the forcing function is 
a gaussian with widths Wx and Wz, respectively, centered at 
{x,z) = {+L/4,0): 



(11) 



Along the inflow (y) direction, the forcing function is antisym- 
metric in y about y — —L/4. It varies linearly with y close to 




the reconnection point and then levels off to a constant value 
(Y — > ±1) over a distance of ^ 2wy upstream of the x-point: 



Y{y) — tanh 



y~L/4 



tanh ( y±I^ 1 + 1 (12) 



The temporal behavior of the forcing is controlled by the 
function Q{t). This function starts at zero and increases 
monotonically with time at the rate dQ/dt ^ 1/t/ = 0.1 
until it plateaus at the value Fao'- 



e{t) = Footanh 



I^ooIt/ 



(13) 



Since the initial slope of Q{t) and hence the initial rate of 
ramping is held constant at 1/r/ 0.1 from one simulation 
to the next, it takes more time to ramp up to the stronger levels 
of forcing (roughly r ^ |-Foo|r/). 

In our simulations the spatial and temporal structure of the 
forcing is controlled by the parameters Wx,Wy,Wz,Tf, and 
Foo. The appropriate choice of these parameters presumably 
varies from system to system, depending on the physical ap- 
plication under consideration and the nature of the forcing 
agent. In this study, as a first step, we examine the recon- 
nection behavior as a function of the asymptotic strength of 
the forcing function Fqo, holding the other forcing function 
parameters fixed at values that are convenient for the size of 
our simulations, e.g. Wx — Wz ~ 5.0di,Tf = 1.25a;~^ . In 
the substantially larger system considered in Ref. |[T9|, the au- 
thors employed a proportionally broader forcing function with 
results similar to those described here. The detailed study of 
how the geometry of the forcing function impacts the recon- 
nection is left for future work. 



III. NUMERICAL RESULTS 

This study includes ten simulations. Five are two dimen- 
sional and differ from each other only in the asymptotic level 
of forcing Foq. These five simulations are included as a base- 
line for comparison with the other five, which are three dimen- 
sional with forcing localized as a gaussian in the out of plane 
direction. In the 3D simulations, the width of the forcing func- 
tion along z is equal to its width along x (wx = Wz — 5.0di). 
In the next section we will describe the results of a typical 2D 
run and then compare to a similar 3D case. First, however, 
some important definitions are required. 

Following Ref. I>14J and our previous 2D work |,19J , we 
focus on the ion inflow Vm (= Vy), outflow Vout (= Vx), 
upstream reconnecting magnetic field Bii (= Bx), and up- 
stream density n^, near the boundaries of the ion dissipation 
region where the ions decouple from both the electrons and 
the magnetic field. As found in earlier work |T9'|, the ion and 
electron velocities typically decouple approximately O.bc/uipi 
upstream of the x-point below the center of the forcing func- 
tion, and we therefore measure the upstream quantities Vin, 
Bd, and rid at this location (+L/4, -L/4 ± 0.5, 0). The re- 
connecting component of the magnetic field does not vary 
strongly as a function of y between 0.5 and 1 .0 ion skin depths 
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upstream of the x-point; therefore, the results presented here 
are not strongly dependent upon the method used to determine 
the upstream edge of the dissipation region. This holds in both 
2D and 3D. 

The electron outflow (i.e. outflow in the x-direction) typ- 
ically becomes quite large very close to the x-point and can 
be in excess of the ion Alfven speed. However, at the down- 
stream edge of the dissipation region, the electrons slow down 
to flow roughly with the ions. The location where the two 
flows come together is typically close to the location of the 
maximum in Vix. Therefore the downstream edge of the dis- 
sipation region is taken to be the location of the maximum in 
Vix- Defining the downstream edge as the location of the ion- 
electron velocity crossover point yields nearly identical results 

miiii. 

A last important definition concerns the reconnection rate. 
In 2D, the reconnection rate is defined as the time derivative 
of the reconnected magnetic flux inside the magnetic island 
per unit length in the (uniform) z direction: 



d 

^o^gll By{x,-Ly/A)dx, 



(14) 



where y = —L/4 is the plane where Bx = 0, and X , X Q are 
the x-locations of the x-point and o-point, respectively. A total 
reconnection rate (rather than a rate per unit length) cannot be 
defined in the 2D limit, since the uniformity of the system 
along z implies that the total reconnected flux is infinite. This 
is not the case in the 3D simulations discussed here, in which 
the extent of the magnetic island and the total reconnected 
flux are both finite in the z direction. Thus in 3D it is natural 
to define the reconnection rate as simply the time derivative of 
the total reconnected flux inside the magnetic island: 



d 

— I Bydxdz 



(15) 



where A is the area of the magnetic island (found numerically 
using field-line tracing) in the central plane of the equilibrium 
current sheet y — —Ly/4. Fig. (jlj) shows a time series By in 
this plane in a simulation discussed in the next section. The 
location of the forcing function, centered at (x, z) = (6.4, 0), 
can be discerned from the dipolar signature in the right half 
of the figures. As we show below, nearly all the reconnected 
flux is inside this roughly circular region under the forcing 
function. (In both 2D and 3D, small secondary islands are 
also typically present; however, the total flux in these islands 
is small and they have no significant impact on the recon- 
nection rates reported here.) The 3D reconnection rate given 



by Eq. (15i has physical units of [flux/time], in contrast to 
Eq. ( [I4| , which has physical units of [flux/time/length]. To 
compare the 2D and 3D reconnection rates, the 3D rate there- 
fore must be somehow translated into a reconnection rate per 
unit length, as in 2D. No unique way exists to do this, mak- 
ing a meaningful comparison of the 2D and 3D reconnection 
rates difficult. Two crude approaches to such a comparison, 
however, are discussed further in the following section. 



A. Temporal Behavior 

Fig. (j2]i shows time series of the dissipation region param- 
eters in typical 2D (left) and 3D (right) simulations. In the 
2D simulation, the asymptotic upstream forcing strength is 
^oo = 0.15, as seen from the solid curve in Fig. (|2^). Note 
the level of forcing has become constant by t « 3. In panel 
(a), we plot the outflow velocity (dashed curve), and 5 times 
the inflow velocity (dotted curve) in the 2D case (the factor of 
5 serves to place the quantities on a similar scale). As the 
forcing is turned on, the inflow velocity increases, causing 
the pressure and density (panel c) to mount under the forc- 
ing function. This back pressure causes the inflow velocity 
to level off and then decrease around t w 4.5. Magnetic flux 
also piles up under the forcing function causing the upstream 
field to increase, as seen in panel (e) during the period t <\2. 
Consistent with past studies {e.g. Refs. |.14, ,16.,191 ), at about 
t ~ 12, when the flux pile-up has narrowed the current sheet 
to about an ion skin depth, the layer opens up into a Petschek- 
like configuration and fast reconnection begins, as can be seen 
by the rise in the reconnection rate i?o in panel (g). This on- 
set of reconnection causes Vin to increase again and B^ to 
level off, indicating that the rate of reconnection is sufficient 
to prevent further pile up of upstream flux, i.e. flux is liberated 
from the dissipation region by reconnection at a rate compa- 
rable to that with which flux is pumped into the region, so 
that Bd is stable for an extended period (t ^ 10 — 25). The 
ratio Vin/Vout ~ 0.1 — 0.2 in this phase is consistent with 
continuity arguments: Assuming in the 2D case that the ion 
diffusion region has a width 25 and a length 2D, the approx- 
imate incompressibility of the ion flow yields VinD ~ Vout^ 
or Vin/Vout ~ 6/ D. As shown in Fig. (j3|l, the value of 5/D 
is typically between 0.1 - 0.2. Here S/D is computed by di- 
viding the half-width at half-max of the current sheet at the 
location of the x-line by the distance from the x-line to the 
location of the peak ion outflow velocity along x. The aspect 
ratio varies somewhat in time, approaching some minimum 
value then generally becoming slightly larger again as the cur- 
rent sheet opens up and shortens into a Petschek-like config- 
uration. Strong reconnection continues in the 2D system for 
about another 10 time units beyond the final time (t = 25) 
shown in the figures, at which point the system runs out of 
magnetic flux to reconnect and the reconnection rate drops 
sharply to zero. The features described here are observed at 
all levels of forcing in 2D that are sufficiently strong to trig- 
ger Alfvenic reconnection. In Ref. |19| qualitatively similar 
results were found in a larger (L x L — 102.4rfi x 102. 4di) 
2D system. Several differences are apparent in comparing the 
2D data shown in left half of Fig. Q with the 3D data in the 
right half. For example, the enhancement in the density (panel 
d) and upstream magnetic field (panel f) are weaker or more 
gradual in 3D than in 2D. This results from a key difference 
between 2D and 3D: in 2D flux can compress upstream of the 
current sheet or it can reconnect, flowing out along x. Simi- 
larly, density can build up under the forcing function or it can 
be cleared out by reconnection. In 3D, plasma and magnetic 
flux have the option of either reconnecting or flowing away 
from the forced zone along the ±z directions. This differ- 
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FIG. 1: (color online) Time evolution of By{x, z)att= 12.1875 (left), t = 16.875 (center), and t = 20.625 (right). The 
reconnected flux moves in the —z direction. Maximum and minimum values of the plotted quantity are shown above each plot. 



ence has at least two important consequences: first, the forc- 
ing function is not as effective at compressing plasma and flux 
in 3D. Indeed, the gradual rise in seen in Fig. (j2]F) results 
largely from the convection of stronger upstream magnetic 
field into the reconnection zone rather than from in-situ com- 
pression. To partially compensate for this, in Fig. ^ we com- 
pare 2D and 3D simulations with somewhat different levels of 
forcing (Foo = 0.2 in 3D vs. F^o = 0.15 in 2D as shown 
in Figs. (j2^,h)). A second consequence, which we discuss 
further below, is that the reconnection "efficiency" in 3D is 
smaller than in 2D - that is, a smaller fraction of the upstream 
magnetic flux that is initially contained within the forced re- 
gion is reconnected in 3D than in 2D, making the reconnecting 
phase of the 3D simulations shorter than those in the 2D runs. 

Two measures of the reconnection rate are shown in 
Fig. (j2ji). The dotted curve labeled Rq, like panel (g), is the 
integral given by Eq. ( 14 1 calculated in the z — plane. This 



quantity is a rough measure of the rate of reconnection im- 
mediately under the center of the forcing function. Although 
one might expect this z = location to yield the largest local 
reconnection rate, in fact substantially larger values can occur 
at negative z values near the edge of the forced zone. This can 
be seen from Fig. Q, which shows the z-dependence of the 
reconnected flux per unit z [Eq. ( 14 1 without the time deriva- 
tive] at various times. The area under the curves at a given 



time yields the total reconnected flux [Eq. (15i without the 
time derivative]. The flux distributions become skewed to the 
left at later times due to the dragging of the magnetic field by 
the electrons. This leftward electron flow carries the bulk of 
the (rightward) electric current density that supports the re- 
versal in the reconnecting magnetic field oc dEx/dy. The 
transport of reconnected flux by the electron flow can also be 
seen in the downward drift of the By profile in the center of 
the sheet shown in Fig. ([T]|. This effect has also been reported 
in other studies of three-dimensional reconnection (see, for 
example, Refs. EOl 1211 l22l |23l). A second measure of the 
reconnection rate is shown as the dashed curve in Fig. (j2}i) la- 
beled i?tot /lO. This curve depicts the total reconnection rate 
Rtot given by Eq. ( 15 i divided by the full width of the forc- 



ing function along z (2wz — lO.Odi). As can be seen from 
Fig. Q, this width roughly characterizes the extent along z of 
the reconnected flux distribution and thus is a crude 

measure of the average reconnection rate per unit length in 
the 3D system. Although the reconnection rates as character- 
ized by these measures are somewhat smaller in the 3D sys- 
tem than those in the 2D case, we show in the the following 
section that, when the 2D and 3D reconnection rates are plot- 
ted as a function of the instantaneous values of the upstream 
reconnecting magnetic field [Figs. (|2^,f)], similar results are 
obtained. That is, for a given value of the upstream magnetic 
field, the rates of reconnection characterized by either Rq or 
Rtot/iO are roughly comparable to the corresponding 2D val- 
ues. Another feature apparent in the 3D data of Fig. Q is the 
downturn in the reconnection rate (panel h) and upstream field 
(panel f) seen att~ 23. At this time, nearly all of the mag- 
netic flux in the forced zone has been expelled by the forcing 
function, and thus the system has effectively run out of flux 
to reconnect. This is much earlier than the otherwise similar 
collapse in the 2D reconnection rate at t ~ 38. As was noted 
earlier, less flux is recoimected in 3D because the plasma and 
magnetic field are free to flow outward along the ±z direc- 
tions without reconnecting. In 2D, the uniformity of the sys- 
tem along z prevents any net transport of flux in that direction, 
and a much larger fraction of the initial upstream flux in 2D 
(over 70% percent compared to about 10% in 3D) ends up 
inside the magnetic islands at late times. 

A final difference between 2D and 3D evident from 
Figs.(|2^,b) is that the ratio Vm/Vout in 3D is enhanced by 
roughly a factor of two over the 2D case. Given that the 
layer aspect ratio 6/D [Fig. (|3j] shows a much smaller in- 
crease, this enhancement is inconsistent with the 2D continu- 
ity argument described earlier, Vin/Vout ^ 5/D. In the 3D 
simulations, however, the continuity relation must be modi- 
fied by a geometric factor In 2D, the ion outflow from the 
reconnection region is almost entirely along the x direction 
with Vix ^ Viz- In 3D, as one can see from the snapshot 
shown in Fig. ([5]l, the ion outflow is nearly omnidirectional 
with Viz ^ Vix, and the ion diffusion region is disk-shaped 
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(like a hockey puck) rather than rectangular. Assuming this 
disk has radius D and half-thickness 6, plasma will thus flow 
into an area Ai^ = ttZ?^ (one surface of the disk) and out 
through an area Aout = 2ttDS. Continuity then demands that 

Vi7iAin — VoutAout or 

= (^) = (?^) ^ (^) 

(16) 

Thus, for the same aspect ratio 6/D, one would expect 



Vin/Vout to be a factor of two larger in the 3D case, as in- 
deed it is. 



B. Scaling of Reconnection 

In this section we test two simple scaling relations for the 
plasma outflow and reconnection rate that have been found in 
past studies to characterize reconnection in some forced and 
unforced 2D systems (e.g. lEl [H [H [HI EH EH [39l HqI). 
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Dissipation Region Aspect Ratio vs. Time 




FIG. 3: Time series of the aspect ratio of the dissipation 
region vs. time in 2D (solid curve) and 3D (dashed curve. 
(Same simulations as the data shown in Fig.(|2])) 
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FIG. 4: Time evolution of the reconnected flux per unit z as a 
function of z (solid curves), and the shape of the forcing 
function(dashed curve). 



In the case of the outflow, the dynamics that cause recon- 
nected field lines to accelerate plasma away from the x-point 
are similar to those of an Alfven wave, and so Vout ^ VAd 
where VAd — Bd/ \/4^7nnin. Regarding the reconnection rate, 
for a given dissipation region of length D and thickness 6, 
the continuity arguments discussed in the last section yield 
V^n ~ Co{S/D)VAd where Co ~ 1 in 2D and Co - 2 in 3D. 
Assuming this inflow along y carries an x-directed magnetic 
field Bd into the dissipation region, the rate of reconnection 
per unit length along z (equivalent to the rate at which mag- 
netic flux enters the reconnection zone per unit length along 



z) is VinBd- We therefore expect: 

B, 



Vout 
i?0 



5_ 
D 

VinBc 



D 



(17a) 
(17b) 



Note the reconnection rate Rq is predicted to scale with the 
square of the upstream field. In normalized units, isolating the 
upstream field Bd as the independent variable, these become 



6_ 
D 

D 



Bd 



1/2 



Bd 



(18a) 
(18b) 



We now test these scaling laws under the assumption that 5/D 
has an approximately constant value (on the order of 0. 1 - 0.2) 
over the duration of the Alfvenic phase of the reconnection 
process. In Figs. (j6^,b), to test the validity of Eq. ( 18a i, the 
quantity Vout \/nd is plotted versus Bd with time as a param- 



eter Each curve represents data from a single simulation and 
is labeled by the value of foo for that simulation. A sixth 2D 
simulation with Foe = 0.025 is included to extend the range 
of data. (This weakest level of forcing is insufficient to pro- 
duce Alfvenic reconnection in 3D, so no 3D analog for this 
simulation is included.) The dotted line of slope unity, plot- 
ted for reference, represents an outflow velocity equal to the 
upstream Alfven speed based on Bd- At the stronger forcing 
levels, the data are indeed in rough agreement with the pre- 
dicted scaling, albeit with significant, order-unity time varia- 
tions within any given simulation. At the weakest forcing lev- 
els, the outflow (like the reconnection rate discussed in a mo- 
ment) falls below the predicted scaling. The forcing strength 
in these simulations is insufficient to narrow the current sheet 

down to the ion skin depth scale, and fast reconnection is 

r~| 1/2 1/4 

never triggered. In Figs. d6t,d), the quantity Rq' nj is plot- 
ted vs. Bd to test the scafing of the reconnection rate given 
by Eq. ( 18b i. Plots based on Rtot/^0 rather than Rq yield 
slightly larger but similar results. The dotted lines represent 
the expected scaling, assuming that Vout/VA = 1-0 and that 
the system has an unvarying dissipation region aspect ratio of 
{5/D) = 1/5 (steeper line) or {5/D) = 1/10 (less steep). 
At the later times when the reconnection rates are highest, the 
data in the more strongly forced simulations roughly follow 
the expected trend. In the most weakly forced runs, however, 
the scaling progressively breaks down; the perturbation of the 
current sheet in these cases is insufficient to trigger fast recon- 
nection before a nearly force-balanced state is reached. 



IV. CONCLUSIONS 

We have examined the scaling behavior of reconnection in 
a forced, three-dimensional, periodic system using two-fluid 
simulations with finite electron inertia. The forcing in the sim- 
ulations was driven by a spatially localized forcing function 
added to the ion momentum equation inside the computational 
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FIG. 5: (color online) The gray scale plot in this figure shows 
Jj; in the central plane of the forced current sheet at 
t = 12.188 in the 3D simulation with = 0.2. The 
over-plotted vectors (blue online) show the relative 
magnitude and direction of = Ji/n in the same plane. 



domain. We initialized the system with a one dimensional, 
broad current sheet equilibrium that was tearing-mode stable; 
three dimensional structure entered through the localization 
of the forcing function. Comparisons were made to analogous 
two-dimensional simulations. 

In the two dimensional case, as found our previous work 
in larger 2D systems |19|, sufficiently strong levels of forc- 
ing were found to produce a quasi-steady Petschek-like re- 
connection configuration with a dissipation region aspect ra- 
tio bjD ^ 0.1 — 0.2. The forcing function produces a pileup 



of magnetic flux in the upstream region and hence an increase 
in the upstream magnetic field, until the reconnection rate be- 
comes sufficient to prevent further pileup of magnetic flux. 
The flux pile-up typically halts once a relatively thin cur- 
rent sheet is formed — ^between 0.5 and 1 .0 ion skin depths in 
width — at which point the magnetic separatrix opens up, and 
fast, Alfvenic reconnection begins, leading to a period when 
the upstream reconnecting field Bd is relatively stable. The 
rate of reconnection per unit length along the (ignorable) z di- 
rection was found to scale roughly like B\, as expected from 
simple scaling arguments, although significant time variations 
are observed during the course of any given run than cannot 
be explained by this simple scaling. 

In three dimensions, freedom of the plasma to flow out in 
the ±z directions makes the forcing function less effective at 
compressing plasma; the resulting flux pile-up is weaker in 
3D, convection dominates over compression in the upstream 
region, and the period of reconnection is shorter than in com- 
parable 2D systems. The localization of the magnetic island 
along z in 3D makes it possible to define a total reconnection 
rate, rather than a rate per unit length, as in 2D. For sufficiently 
strong forcing, this total rate, like in 2D result, was found to 
be roughly proportional B\. 

To directly compare the 2D and 3D reconnection rates one 
must translate the 3D rate into a rate per unit z, as in 2D. Two 
rough methods of extracting such a quasi-2D rate from the 3D 
data were described - one of which is simply to divide the 
total rate in 3D by the width of the forcing function - with 
similar results. When either of these two rates are plotted as a 
function of the upstream magnetic field, the 2D and 3D results 
are comparable. On the other hand, the distribution of recon- 
nected magnetic flux in 3D becomes strongly skewed at late 
times by the flow of electrons in the current layer, and near the 
edge of the forcing function, the flux builds at local rates that 
can substantially exceed the 2D values. 
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